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Abstract. We study integrable models solvable by the nested algebraic Bethe ansatz and 
possessing the GL(3)-invariant i?-matrix. We consider a composite model where the to¬ 
tal monodromy matrix of the model is presented as a product of two partial monodromy 
matrices. Assuming that the last ones can be expanded into series with respect to the in¬ 
verse spectral parameter we calculate matrix elements of the local operators in the basis 
of the transfer matrix eigenstates. We obtain determinant representations for these matrix 
elements. Thus, we solve the inverse scattering problem in a weak sense. 
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1 Introduction 

The algebraic Bethe ansatz was found to be a powerful method for describing the spectrum of 
various quantum integrable models [7, 23, 33, 36]. In this approach, quantum Hamiltonians and 
all other integrals of motion are generated by a transfer matrix. The eigenstates of the latest 
can be found in a systematic way, leading to a set of equations determining the spectrum (Bethe 
equations). 

Despite the significant progress of the algebraic Bethe ansatz in calculating the spectrum, the 
application of this technique to the problem of calculating correlation functions for a long time 
led to much more limited results. It is worth mentioning the papers [10, 11, 21], where series 
representations for correlations of the model of one-dimensional bosons were obtained. Later, 
Eredholm determinant representations for correlation functions of this model were derived by 
the method of dual fields (see [22] and references therein). 

Solution of the quantum inverse scattering problem [20, 25] has opened up new opportunities 
of the algebraic Bethe ansatz. Using this result correlation functions of the XXZ spin chain 
were studied in series of works [8, 9, 15, 18, 19]. The explicit formulas for the local operators 
provided by the inverse scattering problem also had played an important role in calculating their 
form factors [5, 16, 17]. It is worth mentioning, however, that the results of [20, 25] essentially 
were based on the fact that the monodromy matrix of the model could be constructed from the 
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i?-matrix. This is true for various spin chains, but not in general. In the present paper we use 
the approach of [10] for calculating form factors of local operators in quantum GL(3)-invariant 
models. Let us briefly describe the main idea of this method. 

The key equation of the quantum inverse scattering method is the RTT-relation [33, 36] 

Ri2{u,v)Ti{u)T2{v) = T2{v)Ti{u)Ri2{u,v). ( 1 . 1 ) 

Here T{u) is the monodromy matrix, R{u,v) is the R-matrix. In GL(3)-invariant models the 
R-matrix acts in the tensor product of two auxiliary spaces Vi (8) V 2 (14 r\j C^, A; = 1, 2) and has 
the form 

R{u,v) = 1 +g{u,v)P, g{u,v) = —^. (1.2) 

u — V 

Here I is the identity matrix in 14 (8> 14) P is the permutation matrix that exchanges I 4 and 14, 
and c is a constant. The monodromy matrix T{u) acts in (8> R, where R is the Hilbert space 
of the Hamiltonian of the model under consideration. Equation (1.1) holds in the tensor product 
14 (8) 14 ® R) and the matrices Tk{w) act non-trivially in I 4 (8) R. 

The monodromy matrix T{u) of a lattice quantum model is equal to the product of local 
L-operators 

T{u) = Lm{u)---Li{u), (1.3) 

where M is the number of lattice sites, and every R-operator satisfies the RTT-relation with 
the R-matrix (1.2). Gontinuous quantum models appear in the limit M —)• 00 . Let us fix some 
site m (1 < m < M) and define two partial monodromy matrices and as 

r(i)(n) = R^(u)---Li(u), T^^\u) = LM{u)---Lm+i{u). (1.4) 


Then obviously 

r(u) = r( 2 )(u)r(i)(u). (1.5) 

We call such model composite generalized model^ [10]. 

We assume that local L-operators in (1.3) depend on the spectral parameter u as follows^: 

L„(u) = 1-I--L„[0]-I- o(m“^), u—^oo. 
u ^ 

Here 1 is the identity operator in (8> R, and the matrix elements of L„[0] are local operators 
of the model. Then it is easy to see that both partial monodromy matrices T^^\u) have the 
standard expansion over c/u\ 

oo .. 

r«(n) = l +J]r«[n](^)” , Z = l,2, (1.6) 

n=0 

where the partial zero mode r^^^[0] is equal to 

m 

rW[o] = J;l„[o]. (1.7) 

^The authors of [10] used the terminology two-site model. We think that this terminology becomes misleading 
in the case of spin chains. The terminology two-component model used in [34] also becomes misleading in the case 
of multi-component Bose or Fermi gases. We choose to use the terminology introduced in [29]. 

^See [24] for concrete examples of L-operators satisfying the RTT-relation with the GL(3)-invariant R-matrix. 
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In this paper we develop a method of calculating form factors of matrix elements We 

reduce them to the form factors of the monodromy matrix entries Tij(z) studied in our previous 
publications. In this way we obtain determinant representations for the form factors of [0] 
and extract explicitly their dependence on the lattice site number m. Then, taking the lattice 
derivative of the results obtained, we find determinant representations for the form factors of 
the local operators {Lm[0])ij- We do not give here explicit determinant formulas, but the reader 
can find them in [2, 4, 27, 28, 30]. 

The paper is organized as follows. In Section 2 we introduce basic notions of the generalized 
model. In Section 3 we consider composite generalized model. Section 4 contains the main results 
of this paper. There we reduce form factors of local operators to the ones of the monodromy 
matrix entries for the GL(3) case and we conjecture a form for the GL(A^) case. The following 
sections contain the proofs for GL(3). In Section 5 we consider the form factors of the diagonal 
partial zero modes Finally, in Section 6 we study the form factors of the partial zero 

modes for i / j. 


2 Generalized model 


2.1 Bethe vectors 


In the framework of the algebraic Bethe ansatz the entries of T{w) act in a Hilbert space H 
and its dual T-L* that possess a pseudovacuum vector |0) and a dual pseudovacuum vector (0| 
respectively. They are normalized by the condition (OjO) = 1. These vectors are annihilated by 
the operators Tij{w), where i > j for |0) and i < j for (0|. At the same time both vectors are 
eigenvectors for the diagonal entries of the monodromy matrix 


rn('W^)|0) = Ai(u;)|0), (0|rji(r(;) = Ai(t(;)(0|, i = l,2,3, 


where Xi{w) are some scalar functions. In the framework of the generalized model, \i{w) remain 
free functional parameters. Actually, it is always possible to normalize the monodromy matrix 
T{w) —>■ X 2 ^{w)T{w) so as to deal only with the ratios 


ri{w) 


Xijw) 
X2{w) ’ 


rsiw) 


Xajw) 

X2{w)' 


( 2 . 1 ) 


Below we assume that X 2 {w) = 1. 

Bethe vectors are certain polynomials in the operators Tij{u) with i < j acting on the 
pseudovacuum vector |0) [3, 12, 13, 23, 37]. In the GL(3)-invariant models they depend on two 
sets of variables called Bethe parameters. We denote the Bethe vectors Bq ;,(u; u). Here the Bethe 
parameters are u = {ui,..., Ua} and v = {ui,..., Vb}- The subscripts a and b (a, 6 = 0,1,...) 
respectively denote the cardinalities of the sets u and v. 

Similarly we can construct dual Bethe vectors in the dual space as polynomials in the opera¬ 
tors Tij{u) with i> j acting on the dual pseudovacuum vector (Oj. We denote them Ca,b{u\v) 
with the same meaning of the arguments and subscripts. 


2.2 Notation 


Besides the function g{u, v) we also introduce a function /(«, v) 




u — V + c 
u — V 


( 2 . 2 ) 


We denote sets of variables by bar: u, v etc. If necessary, the cardinalities of the sets are 
given in special comments. Individual elements of the sets are denoted by subscripts: Wj, 
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etc. We say that x = x', if #x = ^x' and Xi = x[ (up to a permutation) for i = 1,..., We 
say that x ^ x' otherwise. 

Below we consider partitions of sets into subsets. The notation u => {ni,?Xn} means that the 
set u is divided into two disjoint subsets. As a rule, we use roman numbers for subscripts of 
subsets: Ui, hii etc. However, if we deal with a big quantity of subsets, then we use standard 
arabic numbers for their notation. In such cases we give special comments to avoid ambiguities. 

Similarly to our previous papers (see, for instance, [29]) we use a shorthand notation for 
products of some functions. Namely, if the functions (2.1) or the function / (2.2) depend on 
sets of variables, this means that one should take the product over the corresponding set. For 
example, 

ri{u)= Y{ri{uk), f{z,w)= f{z,Wj), f{u,v)= Yl f{uj,Vk). (2.3) 

Uk£u WjGw Uj£uvk£v 

By definition any product with respect to the empty set is equal to 1. If we have a double 
product, then it is also equal to 1 if at least one of sets is empty. 

In Section 3 we shall introduce several new scalar functions and will extend the conven¬ 
tion (2.3) to their products. 


2.3 On-shell Bethe vectors 


In the algebraic Bethe ansatz the role of a quantum Hamiltonian is played by the transfer matrix. 
It is the trace in the auxiliary space of the monodromy matrix: tr T (u). The eigenstates of the 
transfer matrix are called on-shell Bethe vectors. The eigenstates of the transfer matrix in the 
dual space are called dual on-shell Bethe vectors^. We will denote usual and dual Bethe vectors 
as Ma^b{u,v) and Cafi{u,v) respectively. Different equivalent formulas for these Bethe vectors 
were presented in the paper [3]. In what follows we will not use these explicit expressions. 
Instead, we will use formulas which relate Bethe vectors of the composite model with the ones 
for the components of the model (see formulas (3.4) and (3.6) which are proved in the first 
part of this paper [29]) and the action of the monodromy matrix elements onto Bethe vectors 
obtained in [3]. 

A (dual) Bethe vector becomes on-shell, if the Bethe parameters satisfy the system of Bethe 
equations. We give this system in a slightly unusual form 


ri{ui) 






rsivi) 


/(d,%) 




(2.4) 


These equations should hold for arbitrary partitions of the sets u and v into subsets {-UijtZn} 
and {hi,i;n} respectively. Obviously, it is enough to demand that the system (2.4) is valid for 
the particular case, when the sets Ui and hi consist of only one element. Then it turns into the 
standard system of Bethe equations. 

If the sets u and v satisfy (2.4), then 


tl T{w)B,a,b{u,v) = T{w\u,v)Ma,b{u,v), 


Ca,b{u,v)trT{w) = T{w\u,v)Ca,b{u,v), 


with 


t{w\u, v) = ri{w)f{u, w) + f{w, u)f{v, w) -h r 5 {w)f{w, v). 

Besides usual on-shell Bethe vectors it is also convenient to consider twisted on-shell Bethe 
vectors (see, e.g., [2]). They are eigenstates of a twisted transfer matrix, that in its turn, is the 

®For simplicity here and below we do not distinguish between vectors and dual vectors, because their properties 
are completely analogous to each other. 
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trace in the auxiliary space of the twisted monodromy matrix Tfi{u). The last one is defined as 
Tfi{u) = kT{u), where k = diag(Ki, K 2 , ^ 3 ). The matrix elements Ki {i = 1, 2, 3) are called twist 
parameters. A (dual) Bethe vector becomes twisted on-shell vector, if the Bethe parameters 
satisfy the system of twisted Bethe equations: 


ri{ui) 



fjuuUn) 

f{Uu,Ui) 




rsivi) 



f{Vl,Vu) 




(2.5) 


where ki = and Ui = #Vi. These equations also should hold for arbitrary partitions of the 
sets u and v into subsets. 


2.4 Description of scalar products 

Recall a formula for the scalar product of generic Bethe vectors [32]: 

Sa,b = Ca,fe(h^;h^)Ba,6(M®;u®) = ^ri(Mf)ri(Mj)r3(i;^)r3(uf)/(iif,uJ)/(u®,i2f) (2.6) 






Zau ,61«; < I V ?; Vi ) Zai ,bn «i I <; ) • 


Here all the Bethe parameters are generic complex numbers and the sum is taken over the 
partitions of the sets vP, , and 












{vfX}- 


The partitions are independent except that = Oi with Oi = 0,..., a, and = 

= bi with bi = 0,... ,b. From this we find = On = a — Oi and = 

b^ = b — bi- The functions Za^^bi and .^ai, 6 n so-called highest coefficients. They are equal 
to a partition function of 15-vertex model with special boundary conditions [32]. The reader 
can find their explicit representations in [1, 38]. We do not use these explicit formulas in the 
present paper except Zq^o{0 ] 0 | 0 ; 0 ) = 1. This condition is needed to satisfy the normalization 
5o,o = ( 0 | 0 ) = 1. 

If Cafi{u^',v^) is a twisted on-shell vector and is a usual on-shell vector, then 

we can express the functions ri{vP-'^) and r’ 3 (uj’^) in (2.6) in terms of (twisted) Bethe equa- 

tions (2.4), (2.5). We denote such scalar product by 5^ It is easy to see that 



E 




Kl 


an 


H2\ 






«3/ 


X f{v^,U^)f{vf,uf)Zaj, (uJ; Un I Vi ; vf ) Za, ,6jj (uf; <|). 


(2.7) 


Remark 2.1. One should be careful using Bethe equations in the sums over partitions, be¬ 
cause there might be problems if some parameters from the sets u'^, coincide with the ones 
from M®, (i.e., 7 ^ 0 or/and v^^Civ^ 7 ^ 0 ). The matter is that after imposing the Bethe 

equations we cannot consider the limit where one solution of Bethe equations goes to another 
solution. Instead one should first take the limit in (2.6) (what leads to the appearance of the 
derivatives r[{u) and P^{v)) and only then impose Bethe equations. However, if we consider 
the scalar product of twisted and usual on-shell vectors, then we can use the twisted and usual 
Bethe equations from the very beginning. The matter is that in this case the parameters 
and are functions of k = { 1 ^ 1 , K 2 , K 3 }: = u^{k) and = v^{k). Therefore we always 

can take k such that vP = 0 and H h® = 0 . Then, if necessary, we can consider the 

limit where some parameters coincide. In this case we should treat the variables vP and as 
functions of k. 
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The sum over partitions (2.7) was studied in [2, 4] for arbitrary \dh\QS of the Bethe parameters. 
There this sum was reduced to a single determinant in the case ki = K 3 . It was proved that the 
determinant vanishes at"^ R = 1. Thus, if 7 ^ 0 or 7 ^ 0 , then setting it = 1 in (2.7) we 
obtain an identity 


0 = 


^ (^n ) l^i i )'^ai,6n (^i ) l^n ) )• 


( 2 . 8 ) 


In the exceptional case = 0 we have 



= 1 . 

/t=i 


(2.9) 


Remark 2.2. One should not be surprised that the r.h.s. of (2.8) does not give the norm of an 
on-shell Bethe vector in the case {u^,v^} = {u^,v^}. Indeed, as we explained above, in order 
to obtain the norm one should consider the limit {u^(R), (R)} —>■ {u^,v^} at it —>■ 1 in (2.7). 
Instead we simply set k = 1 in (2.7) for generic values of the Bethe parameters. In this case 
we obtain that the sum ( 2 . 8 ) vanishes [ 2 ], and then this result can be continued to the point 

Identity (2.8) plays a central role in the calculation of form factors of local operators. 


2.5 Universal form factors 

Form factors of the monodromy matrix entries are defined as 

Tifiz) = = C“'’''(u^;i;^)U,-(z)B“’'(u^;i;^), 

where both C°''{vF]v^) and B“’^('u®; u^) are on-shell Bethe vectors, and 
a' = a + 5ii — 5ji, b' = b + 5j^ — < 5 * 3 . 


The parameter 2 ; is an arbitrary complex number. 

It was proved in [30] that if {u^,v^} 7 ^ {u^,v^}, then the combination 




v^; M®,u®) 
t{z I u^,v^) 


does not depend on 2 . We call {u^,v^;u^,v^) the universal form factor of the operator 

Tij(z). li Ciu^ = 0 and Civ^ = 0 , then the universal form factor is determined by the 
i?-matrix only. It does not depend on a specific model, in particular, on the functions ri( 2 ;) 
and r 3 (z). 


3 Composite generalized model 

Consider a composite generalized model defined by (1.4), (1.5). Every TF(^u) satisfies RTT- 
relation (1.1) and has its own vacuum state Hereby |0) = ® The opera¬ 

tors Tp\u) and Tp\v) commute with each other, as they act in different spaces. 

Let 

7;«(n)|0)« = Af^(u)|0)«, / = 1,2. 

■^Here and below the notation ii — 1 means Ki = 1, i = 1, 2, 3. 


(3.1) 
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We also introduce 


rl>{u) = 


\f{u) 


1 = 1,2, k = l,3. 


Obviously 

Xi{u) = X[^\u)xf\u), rk{u) = r^^\u)r^^\u). 

Below we express form factors in terms of (u), therefore we introduce a special notation for 
these functions 

rk{u) 


r^^\u) = £k{u), and hence, 




k = l,3. 


Observe that (1.6) implies 
c 

■ u 

and due to (1.5) 

u 


4(m) = 1+ 4[o]-+ o(tt ^), 


rk{u) = l + rfc[0]- + o(u ^). 


(3.2) 

(3.3) 


We extend convention (2.3) to the products of the functions r^\u) and ^fe(tt). Namely, whenever 
these functions depend on sets of variables (for instance r^l\u) or ik{v-a)) this means the product 
over the corresponding set. 

3.1 Bethe vectors and partial Bethe vectors 

We can introduce partial Bethe vectors B^^]^(u;'i;) for both partial monodromy matrices T^^\u). 
Then Bethe vectors of the total monodromy matrix can be expressed in terms of the partial 
Bethe vectors as follows [6, 29]: 




(3.4) 


The sum is taken over partitions u => {ui, ttn} and v => {ui, Un}. The cardinalities of the subsets 
are given by the subscripts of the Bethe vectors. 

(2i^ 1 

We can present the product of functions r) ^4) as r\ (uj) = ri(4)4 see (3.2). More¬ 
over, if we deal with an on-shell Bethe vector, we can express ri('Ui) in terms of the function /, 
thanks to the Bethe equations (2.4). Then we obtain 


Ma,b{u-,v) = ^ ^^||5y/(t6i,Wn)/(Un,Ui)/(Ui,Ui)B^^^^^(ni;Ui)B^^^^_j,^(Un;Un)- 

Similarly, dual Bethe vectors can be expressed in terms of partial dual Bethe vectors 


(3.5) 


(3.6) 

where the sum is taken again over partitions u {tii, tin} and v {4,4}- 

If Ca,b(u;v) is a twisted on-shell Bethe vector, then we can use again (3.2) as 44 ) = 
r3(4)4^(4) and express ^ 3 ( 4 ) through the twisted Bethe equations (2.5). We get 


^15(4^) = E (“) ^^/(^i>^n)/(?^n,4)/(^^i,4)CiE(4Wi)Ci^E(^“’^“)- 

Here we have added the superscript (k) to the vector C^^^^(n;i;) in order to stress that it is 
a twisted dual on-shell Bethe vector. 
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3.2 The action of total and partial zero modes 

The action of the operators Tj^j\z) on the corresponding partial Bethe vectors is the 

same as the action of total Tij{z) on the total Bethe vectors v) [3]. One should only replace 

in the formulas the functions rk{z) by their partial analogs r® (z). The same replacement should 
be done in the action of the partial zero modes on the partial Bethe vectors [30] . In this section 
we give some of those actions used below. 

The action of the total zero modes Tij[0] (with i < j) on the total Bethe vectors can be 
easily extracted from the formulas given in Appendix A of the first part of this paper [29] using 
expansion of the monodromy matrix elements (1.6). They are 


Ti 3 [ 0 ]Ba b(u;u) = lim ^Ma+i,b+i{{'w,u}-, {w,v}), 

W^OO ^ 

Ti 2 [ 0 ]Ba b(u;i;) = lim ^Ba+i,fe({rc, u}; u), 
T 23 [ 0 ]Ba,b(tt;i;) = lim ^Ba,fe+i(-u; {tc, u}). 


(3.8) 


The right action of the operators T)i[0] with i < j on dual Bethe vectors is quite analogous. One 
should replace in (3.8) ^^[O] by Tji[Q] and Mafi{u]v) by Ca^biu]v). 

The action of the partial zero modes Tjj^^[0] (with i < j) on the partial Bethe vectors B^^^ is 
similar to (3.8): 

^1? [0]®S(^; v) = 1^ f Bi^i,b+i({w^, u}] {u;, u}), 


Tu[opaii^;v) = 1 ^ 


W^OO 

The action of the partial zero modes [0] has the following form 
= (4[0] - a)B2(«;i;), 

^ 2 ?[o]®S(^;^) = (« - 

+ (3.10) 

where ik[0] are determined by (3.3). In all the formulas above Bethe vectors (partial or total) 
are generic. 

In Section 6 we will also use singular properties of on-shell (dual) Bethe vectors 

Ca,b{u^-,v^)Tij[0] = 0, Tji[0]Ea,b{u^',v^) = 0, i < j. (3.11) 

Here Ca,b{u ^and Ba,fe(tt‘^; u'^) are on-shell Bethe vectors. This property was found in [26] 
for GL(A^)-invariant models. In the GL(3) case it also follows from the explicit formulas of the 
action of the operators Tij{z) onto Bethe vectors [3]. 


(3.9) 


4 Main results 


Theorem 4.1. Let Ca\b'{u^]v^) andMafi{u^;v^) be total on-shell vectors such that {u^^ 
{that is, these on-shell vectors have different eigenvalues). Then 


Ca\b’iu^ 



\il{u^)£3iv^) 





(4.1) 


where is the universal form factor of the total operator Tij{z) and a' = a Sn — 6ji, 

b' = b-\-6j3 - Si3. 
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Theorem 4.2. Let Ma^b{u',v) be a total on-shell vector and Ca,b{u',v) its dual on-shell vector. 
Let Ca,b{u{R);v{R)) he a deformation ofCafiiu',v) such that the parameters u{R) andv{R) satisfy 
twisted Bethe equations (2.5), and u{R) = u, v{R) = v at R = 1. Then 


:a,b{u]v)T^P 


,^b{u] v) = + <5i,3£3[0] + ^ log 

\ dKi Is[v{k)) 


K =1 J 


ia,b 


u:v 


\i4.2) 


The proofs of these theorems will be given in the next sections. 

If the partial monodromy matrix has the structure (1.4), then the functions ik{u) 

actually depend also on the number m: 


m 

ik{u) = Ylhiuln), k = l,3, 

n=l 


where ik{u\n) are the local ratios 
Ai(ri|n) 


ii{u\n) = 


Q(ri|n) = 


hiu\n) 


A2(ri|n)’ A2(n|n) 

In (4.3), we introduced the vacuum eigenvalues of local L-operators Ln{u) 


(4.3) 


{Ln{u))ii\0) = Aj(tt|n)|0), i = 1, 2, 3. 


Using (1.7) and Theorems 4.1, 4.2 we can find form factors of the local operators {Lm[0])ij, 
i,j = 1,2,3. Namely, one has simply to consider the difference of two T^^\u) based on m 
and m — 1 respectively. 

If {u^, v^} / {-u®, u^}, then we have 


Ca',b' (^i^; t’^) (Lm[0]) .^.]Ba_b(n®; u®) 

\ii{u^\m)l‘s{v^\m) ) y^i(u^|n)73(i;'^|n) j 





If {n*^, v^} = {u®, v^} = {n, u}, then 


Ca,b('U; v) {Lmff\) ^Ma,b{u] u) 


d ^ ii[u{K)\m) 
dKi l^{y{R)\m') 


K=1 


2,b{u;v 




If we deal with a continuum model, then form factors of local operators can be found directly 
from (4.1), (4.2). In this case the integer number m turns into a continuous variable x. This 
parameter enters only the functions ik, and taking the x-derivative of (4.1), (4.2) we find form 
factors of local operators in the point x. 

Thus, we obtain form factors of local operators in the generalized model without use of 
a specific representation of the algebra (1.1). In fact, it means that we have a solution of the 
quantum inverse scattering problem in the weak sense. We cannot express the local operators 
in terms of the monodromy matrix entries as it was done in [20, 25], but we can find all 
their matrix elements in the basis of the transfer matrix eigenstates. Furthermore, we have 
determinant formulas for all these matrix elements [2, 4, 27, 28, 30]. 

GL(A/') generalisation. We would like to mention that Theorems 4.1 and 4.2 admit a direct 
generalization to GL(A^)-invariant models with N > 3. Indeed, Bethe vectors (and dual ones) 
of GL(iV) models depend on N" — 1 sets of parameters U = {fj, • ■ •, j = Ij 2,..., N" — 1 
and N — 1 integers aj that correspond to the cardinalities of each set: 

]Ba(t) = IBai,a2,...,ajv-i > Ca(t) = Cai,a2,...,ajv-i • 
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The action of the diagonal entries on the vacuum vectors is similar to (3.1) 


TP\t)\0f^ = xf\t)\0f\ 1 = 1,2, i = l,...,iV, 


and we can introduce 


ai{t) 


Ai‘'W 

A'hW’ 


i = l,...,N -1. 


Note that in the case N = 3 we have ai{t) = £i{t), while a 2 {t) = ^(t). 

Conjecture 4.3. Form factors of the partial zero modes T^p[0] in GL(A^)-in?;ariant models are 
given by 


C-a{t)T^Pm-a{t) 




n 1 j ^ ’ for i 


where is the universal {z-independent) form faetor of the total operator Tij{z) and we 

extended the convention (2.3) to the funetions oik- t{K) is a set of n-twisted on-shell Bethe 
parameters, coinciding with t when R = 1. 


This conjecture generalizes Theorems 4.1 and 4.2, proved for N = 2 and N = 3. 

Remark that since C5(s)Tjj[0]]Ba(t) = 0 when s ^ t, the above conjecture and theorems 
provide also the form factors for Tk [0]. 

It is worth mentioning that in the cases N = 2 and N = 3, compact determinant representa¬ 
tions for the universal form factors are known. In contrast, in the case N > 3 such representations 
are missing up to now. Nevertheless, if Conjecture 4.3 is valid in the GL(A^) case, then it gives 
explicit dependence on the lattice site m of the partial zero modes form factors. 


5 Form factors of diagonal operators 


We begin our consideration with the form factors of the diagonal partial zero modes t//^[ 0]. It 
is convenient to construct a special generating functional for these form factors [10]. Consider 
an operator 


2=1 


where /3j are some complex numbers. The generating functional is 


(5.1) 


Here lBa,b{u^]v^) is an on-shell Bethe vector, {u^; v^) is a dual twisted on-shell Bethe vector 
with the twist parameters Ki = . 

Lemma 5.1. Let is defined as in (5.1). Then 


a,b a,b^ 

iR'\ 

where is the scalar produet of the twisted and the usual on-shell Bethe veetors (2.7). 


(5.2) 
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It is worth mentioning that an analog of (5.2) for GL(2)-based models was obtained in [14]. 
We will give a proof of Lemma 5.1 in Section 5.1. Now we show how equation (5.2) implies some 
statements of Theorems 4.1 and 4.2. 

Differentiating (5.1) over k, at = 1 we obtain 



Pay attention that the dual vector Ca,b{u^'iV^) = u*^) 

the second term in the r.h.s. of (5.3) is a form factor of the partial zero mode rP''[0]. 
On the other hand, differentiating the r.h.s. of (5.2) over Kj at it = 1 we find 


K,= l 


is an on-shell vector. Thus, 
(Pr 


[o]+/33^3[ 0 ]c(«) 
dKi ii{u^)£3{v^) 



d ^( k ) 

K=1 ii{u^)£3{v^) 

K=idKi 


+ ' 


5MQ[0]+5.,34[0] + /-log^ 
dKi h{v) 


K=1 


2,b{'^;v) 


(5.4) 


where 


1 , if = 

0 , if {u^,v^}. 


Deriving this formula we used the orthogonality of on-shell Bethe vectors depending on different 
Bethe parameters. Comparing equations (5.3) and (5.4) we immediately arrive at the statement 
of Theorem 4.2. If {u'^, ^ {u^,v^}, then we obtain 


,b{u^-v^)T^P[0]Ma,b{u^-,v^) = 


ll{u^)i3{v'~ 


- 1 


R=i dm 


n=l 


It was proved in [4] that the Kj-derivative of the scalar product at k = 1 is equal to the 
universal form factor of the operator Ta^z). Thus, the statement of Theorem 4.1 is proved for 
partial zero modes t//^[0 ]. 


5.1 Proof of Lemma 5.1 


The proof of Lemma 5.1 is lengthy but straightforward. Let us first sketch the general strategy 
before going into details. Knowing the action of the operator on the partial Bethe vectors we 
find its action on the total Bethe vectors. Then we can can calculate the matrix element M 'd in 

a,b 

terms of scalar products of partial Bethe vectors, for which we use equation (2.6). The resulting 
formula becomes rather cumbersome. In particular, it contains a sum over partitions of every 
set of the original Bethe parameters into four subsets. Therefore in this section we use standard 
arabic indices in order to label these subsets. New subsets of the Bethe parameters can be easily 
recombined into new sets of variables, and after this the proof reduces to the use of identity (2.8) 
and equation (2.7). 

Let us now give the details. Using (3.10) and (3.5) we find the action of on the total 
on-shell Bethe vector 




E 


= > 6 ' 


o/3i (^1 [0] - ai) -1-/32 (ai - fei) -1-/33 (^3 [0]+6i) 


^3«) 


X /«, 0 /«wf)/«,< 


Q(nf) 


an, On 
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Multiplying this equation from the left by and using (3.7) we obtain 


nfi^) ^ V ^/3i(€ifO]-aT)+feaT+/33fefO] ^l(^n )^3«) 

^ h{uf)£3{v^) 

X f{uf,u^)f{uf,Un)f{vn, Vi)f{vE,vnf{vi, u?) 






(5.5) 


(r) 

Thus, we have obtained the expression for ^ in terms of scalar products of partial Bethe 
vectors. Note that in spite of the sets and 7® satisfy the Bethe equations (2.4), and the 
sets and v'^ satisfy the twisted Bethe equations (2.5), the partial Bethe vectors in (5.5) are 
not (twisted) on-shell vectors. In other words we deal with the scalar products of generic Bethe 
vectors in (5.5). Therefore, we do not write the additional superscript (k) for the dual vectors 
and we should use (2.6) for the calculation of their scalar products. Hereby, for the scalar product 
of the vectors and we should replace in (2.6) the functions by Ik, while for the scalar 

product of the vectors and B^^) we should replace in (2.6) the functions ru by rk£^^■ 

The use of (2.6) introduces new partitions of the subsets of Bethe parameters, so that, as 
mentioned above (see Section 2.2), we use now arabic numbers to label the numerous subsubsets. 
Thus, we have for the first scalar product 


C^;^«;i;f)B^;^(nf;i;f) = ^£i(kf)^i(kf)4(i;3^)^3(ilf)/(«f,«?)/(nf,ll 


kf) 


X f{v3, )/(?^f, (^3 ; «3 kf; )^ai,fe3 «; Ui\v^ ] U3 ). 

/m jj\yi ) “i ) 

The summation is taken with respect to the partitions 


ttr 


c-C,B -C,B'i 

{''^T }> 




r-C,B -C.B'i 

,^3 I 


The cardinalities of the subsubsets are bn = n = 1, 3. 

Similarly 


n(wf)n«)^’3(^'4)^3(^'f) 




/(M^,Mj)/«,uf; 


(5.6) 


X /«, )/(uf, uf ) ^a4 ,b2iu4-,uf\v^;v^)Za2,bA^2] V^] 


Here the sum is taken over partitions 


ur 


{u. 


-C,B -C,B 


U. 


% 


yC,B 


{v. 


C,B -C,B 


}• 


The cardinalities of the subsubsets are still denoted by On = #Un'^, and bn = i^Vn'^, n = 2,4. 

Now we should express the products of the functions in (5.6) via the (twisted) Bethe 
equations for the full sets {u^,v^} and {u^,v^}. We have 


/(u® 


2 ) 


fiuf, 

fivf, 


n{u2) 
rsiv^) 

ri(kj) 

r^iv^) = 


Z^J[V ,U2), 


U2)fiui,U^)f{uf,ui 


= e 


f{v, 

ai{P2 


f{v2,U^), 


0/(^^f,'*^f)/(^'f,0' 

-0,) fiui,uf)fiu^,U^)f{uf,U^) , c -CN 

f{u^4^ul)f{uq,U%)f{ui,U<iy^ ’ 

/33) /K^^4)/(^2>^4)/(^3^^4) f(-C -CN 

m:vf)f{vf,v^)f{vyv^y^ 
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All these expressions should be substituted into (5.5). After simple but exhausting algebra 
we obtain 


^ ft ^1 f0l+/3.sfa f0]+(fe-/3i )(a-aA+(fe-feife4 (^2 )^1 (^3)^3 (^’f )Q «) 

^ Q(hf)£i(hf)Q(p)£3(i;f) 


X Z 

^ UU-^ VV^ VU-^ UU-^ VV^ VU^ ' 


(5.7) 


Here the sum is taken over partitions of every set of Bethe parameters into four subsets 




t-C,B -C.B -U.B -U.B'i -C.B , f-C.B -U.B -U.B -C.B'i 

{'^1 ) '^2 j ^3 ’^4 }> ^ >'^2 ’^3 ’’^4 }■ 


We have #ix® = and = bn, n = 1,... ,4, but the values a„ and bn 

are free. Note that 02 and 64 explicitly appear as coefficients in (5.7), so that manipulations 
with U 2 ’^ and should include these coefficients. 

The factor Z in (5.7) is the product of four highest coefficients 

^ ,61 (hs ; hf I hf; hf) ,b3 (iif; uf I hf; hf) ,62 (mJ; uf I ; uf) ,64 (uf I uf; ). 

The other factors in (5.7) denoted by F with different subscripts and superscripts are products 
of / functions: 

pZ = fiu2,ui)f{u^,u2)f{u'i,u^)f{u‘i,u^)f{u'{,u^)f{u^,u^), 
pZ = /(hf,hf)/(hf,hf)/(hf,Mf)/(hf,hf)/(hf,?2f)/(uf,hf), 

Pvv = fiv?,vT)f{v2,V3)f{v^,V^)f{v^,V^)f{v^,V^)f{v^,Vi), 

Pvv = f{vt,vf)f{v^,v 2 )f{vf,v^)f{vl,v^)f{vf,v^)f{vl,v^), 

Pyu = fiv?,u'i)f{v^,u'i)f{v^,u'{)fiv^,u'{)f{v^,u'i)f{v^,u^), 

Pyu = f{v3,ui)f{v^,U^)f{vi,U^)f{v^,ui)f{vf,U^)f{v^,uf). 

It remains to combine the subsubsets into new groups: 


;r.C,B 


f-U,B -U.B'I -C 

{uL ,uF I = u- 

r-C,B -C.B'I -C.B 
jUl’ ,U4’ } = Ui ’ , 


r-C,B -C.B'I -C.B 
1% }=^ii > 

r-C,B -C.B'I -C,B 

W2 W3 I = w • 


Then we recast (5.7) as follows: 

^l(«i4)^3(hf 


= 

a,b 


E 


-f{uf,ufi)f (ufi, uf )/(?;?, uf)/(uf, Ujf ^ 


uCiB^{nC,B.C,B^ 

CiB^{n'^iB-C,B^ 




X fivf, uf)f{vy,,Uy,)Gi{uf, uf; u?, Ujf )G2«, <; uf, vf), 
where factors Gi and G 2 are given as sums over partitions 


(5.8) 


Gi«,uf;Ui'?,Uif) = 


E 


/«,<)/«, 0/(^'2> ^^3 )/(^3 Wf 


-C,B^ (-C,B -C,B^ 
-C,B^ f-C,B -C,B 


} 


U2 ,t/3 

X / (113 > ^4 ) / (^2 : ) ■2'a4 ,b2 (^4 wf I ^2 ; ^2 ) ^ai ,63 


1! Ts; '*^3 ^ 


(5.9) 


G2('u?,u(?;uf,uf) = 


E 


o0ih[O]-\-133^3 [0]+(^2-/5i)(a-a2)+(^2-/33)^4 


-C,B^ (-C,B -C,B~, 
=^{^2 ’^3 J 

-C,B^ f-C,B -C,B~, 
U. ,114 ’ } 


and 
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X fiu3,U2)fiu2,ui)f{vf, V^)fiv^,vf)f{v2,U^)f{vf,U^) 

X Za3^bi{u3;ui\vf;vf)Za2,bi{u2',U2\v^;v^)- (5.10) 

It is not difficult to see that the sum over partitions in (5.9) coincides with the sum in (2.8) up 
to relabeling of the subsets: —>■ uf’^, —>■ and —>■ ■ Thus, 

we conclude that Gi = 0 unless uf'^ = 0 and = 0 . Hence, ■u?’'® = hf’'® = v^’^, and 

ai = 04 = 0, 62 = 53 = 0. Then due to (2.9) Gi = 1. 

Looking now at (5.10) and comparing it with (2.7) we see that they coincide up to the common 
prefactor relabeling of the subsets: uf’^, —>• hn’®, —>• 

and —>■ hj’®. Hence, 

_ g^Rl[0]+/334[0]^(«)_ 

Substituting this into (5.8) and setting there uf'^ = 0 and h?’® = 0 we immediately arrive 
at (5.2). 


6 Form factor for off-diagonal partial zero modes 


Now we study the form factors of off-diagonal partial zero modes. We apply a strategy similar 
to the one used in [27], using the commutation relations of the zero modes and properties of the 
type (3.9), (3.10). 

First, we note that the i?TT-relation (1.1) implies in particular for the (partial) zero modes 

[T,,[o],r,-,[o]] =r,-,[o]. 

, 0 . 

II 


[Tij [0], Tii [Oj] = Tij [0], 

, 0 . 

II 


[Tij [ 0 ], Tfcj [ 0 ]] = Tkj [ 0 ], 

[TP[0],TPm = Tp[0]. 

, 1 = 1,2. 

Now, using ni[0] = tP\0]+tP'>[0] and [i;^^[0],r(f)[0]] 

= 0 we conclude 

[ 4 ')[o],r,,[o]] = r;P[o], 



mj[0lTp\0]] = Tp\0], 

i + h 


m0lTp\0]] = 4f[0], 


( 6 . 1 ) 


which are the central relations that we will use for our calculations. 
As a notation, we will note the form factor of Tj^p [0] as 



with a' = a + 6 ii — 6 ji and b' = b + Sj3 — 6 i 3 , and where both vectors are on-shell. 


6.1 Form factor of Ti2^[0] 

We start with Mpp\ the form factor of Tg 2 ^[ 0 ]. We make the calculation in details, the other 

( 12 ) (2 2 ) 

ones following the same steps. To get ^ , we start with the form factor and send one 

of the parameters in to infinity (keeping all the other parameters finite), using relations (3.8): 


w 


w^oo c 


lim -Ca+i,fe(tf^; [ 0 ]Ba+i,fe({rt®, u;}; v^) 

w^oo c 

Ca+i,fe(i2^;h^)rW[0]ri2[0]B„,fe(u^;i;^). 
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Now, since Ca+i^biu^]v^) is on-shell with all parameters finite, it is annihilated by ri2[0] due 
to (3.11). Thus, we can replace the product T22^[0]Ti2[0] by its commutator which in turn gives 
T^^2^[0] through (6.1). It leads to 


w^oo c 


( 6 . 2 ) 


( 2 , 2 ) 

It remains to compute the limit of the explicit expression (4.1) for that was proved in 

the previous section. It is obvious that lim ikiw) = 1, = 1,3, and it has been shown in [27] 

W^OO 

that^ 


lim ,v^) 

w^oo c 


(6.3) 


Hence, we get expression (4.1) for 


(L2) 


6.2 Other form factors 

As already mentioned, the calculation for other form factors follows the same steps, so we just 
sketch the proofs. 

Form factor of T23^[0]. To get we start with ^w}) and take 

the limit w —)■ oo. It makes appear T23[0] that annihilates Ca,fe+i('it‘^;u*^). It leads to 

lim = (6.4) 

w^oo C 

Form factor of T2i^[0]. To get we start with and take 

the limit w —)• oo. It makes appear r2i[0] on the left, and it annihilates ',v^) due to (3.11). 

Then, we can again replace the product T2i[0]T22^[0] by its commutator, and we obtain 


lim ,v^-u^,v^). 

w^oo c ’ ’ 

Form factor of T’32'^[0]. The form factor ' is obtained through the limit 

w^oo c ’ ’ 


(6.5) 


( 6 . 6 ) 


Form factor 


of T^^ 3 ^ [0]. To get , it is convenient to start with the already known form 

(2 3 ) 

factor ^w},v^). Taking the limit tc —>• oo we obtain ri2[0] that annihilates 

Ca+i^;,(n‘^; h'^). Hence, due to the last equation (6.1) we obtain 


r(L3) • 


W—^(X) c u.,o 


- /\/r(h3)/-c ;r,c.r,i 


(6.7) 


Form factor of T3^^[0]. To get we start with {{u'^,v^) and take 


the limit w ^ oo. We obtain 

C ’ ’ 


( 6 . 8 ) 


(2 2) (i j) 

Thus, starting from one initial form factor ^ ' we can obtain all other form factors 

in the special limits of the Bethe parameters. This property is a direct consequence of the 

property of the form factors of the monodromy matrix entries [30]. 


^Strictly speaking the proof of [28] was done for the full form factor, bnt it extends straightforwardly to the 
universal form factor. 
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Consistency with morphisms. It is worth mentioning that there exist also other relations 
between different form factors [28, 30]. These relations appear due to morphisms of the 

algebra (1.1). It was shown in [3] that the mappings 


ip: Tij{u) ^ Tji{u), ip: Tij{u) ^ T 4 _j^ 4 _i{-u), 


(6.9) 


define morphisms of the algebra (1.1). Hereby, the mapping y? is an isomorphism, while '0 is an 
antimorphism. Both mappings (6.9) exchange the partial monodromy matrices of the composite 
model [29] 

{u) ^ Tt^Ui-n)■ 

Transformations (6.9) induce relations between different form factors of the monodromy 
matrix entries (see [27, 28] for details). The latest, in their turn, give us relations between 

U j) 

the universal form factors 


-c», -s»). 


Thus, actually it is enough to compute only four form factors of the partial zero modes 


( 6 . 10 ) 

TifiO] 


(j = 1,2,3) and r22^[0]. All other form factors can be obtained by the mappings described 
above. Using the explicit determinant representations for the from factors one can check that 
they enjoy both the limiting procedures (6.2), (6.4)-(6.8), and the transformations (6.10). 

Generalisation to GL(Al). One can apply the same procedure to the GL(A^) case. In 
particular, the commutation relations (6.1) are still valid, and the limits of the type (3.9), (3.10) 
and (6.3) have been proven for GL(A’) in [27]. The singular vector properties of the type (3.11) 
were obtained for Gh{N) in [26]. All that allows us to relate the different form factors in the 
same way we did for GL(3). It is easy to check that the Conjecture (4.3) is consistent with these 
relations. 


7 Conclusion 

In this paper we have studied form factors of the partial zero modes in a composite generalized 
model with GL(3)-invariant R-matrix. We have reduced these form factors to the ones of the 
monodromy matrix entries Tij{u) considered in our previous publications. As we have mentioned 
already, it means that we have a solution of the inverse scattering problem in the weak sense. 
Apparently the same type of the formulas remains true for GL(A^)-invariant models with N > 3. 

We are planning to apply these results to the study of form factors of local operators in the 
model of one-dimensional two-component Bose gas with d-function interaction [35, 39]. Several 
form factors in this model were studied already in the framework of the coordinate Bethe ansatz 
for some particular cases of Bethe vectors [31]. We are going to apply the algebraic Bethe ansatz 
techniques in order to compute the form factors in the general case. This model possesses GL(3)- 
invariant R-matrix (1.2) and can be considered in the framework of the scheme described in the 
present paper. However the asymptotic expansion of the monodromy matrices (1.6) in the case 
of Bose gas should be modified. In its turn, this modification leads to a modification of the 
partial zero modes. Nevertheless the derivation of determinant representations for the form 
factors of the diagonal partial zero modes t//^[ 0] does not change. A possibility to use this 
result for obtaining all other form factors is not obvious, however they can be calculated by 
a straightforward method similar to the one that we used for the calculation of t//^[0 ] in the 
present paper. This will be the subject of our further publication. 
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